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Abstract 

The algebraic area probability distribution of closed planar random walks of length on a 
square lattice is considered. The generating function for the distribution satisfies a recurrence 
relation in which the combinatorics is encoded. A particular case generalizes the g-binomial 
theorem to the case of three addends. The distribution fits the Levy probability distribution for 
Brownian curves with its first-order 1/A^ correction quite well, even for A^ rather small. 
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1 Introduction 

We are considering closed planar random walks on a square lattice. We aim at finding the 
probability distribution of the algebraic area A enclosed by a random walk of length A^, starting 
and ending at the origin. 

Random walks are a discretized version of closed continuous Brownian curves, which are 
obtained in the limit N ^ oo. In this limit the average size of the walk diverges as VN, its 
average area is proportional to A^, therefore its finite renormalized area is a = A/N. As first 
shown by Levy pj, the asymptotic probability distribution is 

vr 



cosh (27ra) 

Discrete random walks whose links can point in arbitrary directions were previously considered 
[3, and the Levy distribution recovered in the continuous limit [3J. Certain analytic results were 
also obtained concerning the area distribution of walks in the presence of random traps [1], of 
directed random walks (those that begin and end on the y = line) [5], as well as walks confined 
to a finite-size box [6]. 
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The problem of the random walks area distribution and n-winding sectors area distribution 
on a square lattice arises, for example, in the context of random magnetic impurities and the 
integer quantum Hall effect [7]. Also, the distribution is connected to the partition function of 
a lattice gauge model with Z(2) gauge group interacting with a Z(2)-valued Higgs field [8]. 

The problem has been attacked from two sides. For big A^, a finite-size correction to the 
Levy distribution at first order in was derived in Ref. [9] [see Eq. (f29ll below] by relating 
the number of walks to the trace of the Hamiltonian of the Harper model [10] . which is in turn 
related to the Hofstadter model [H] of an electron moving on a two-dimensional square lattice 
in the presence of a uniform magnetic field orthogonal to the plane. On the other hand, for a 
finite A^, all the probabilities involved are rational numbers, therefore the sought distribution is 
a rational function. Properties of its k-th moment — which turns out to be a rational function 
of A^ with integer coefficients — have been studied in Ref. [I2j, where the Levy distribution was 
also explicitly obtained in the A^ — > oo limit (see also [13] )• An exact expression for the finite A^ 
distribution has, however, remained out of reach. 

In this paper we derive a recurrence relation for the generating function of the area probability 
distribution for an arbitrary A^, which is interpreted within a statistical mechanics approach as 
well as in terms of g-commuting operators. The asymptotic limit of the distribution and its 
first-order 1/A^-correction are reproduced numerically. 

2 Generating function of the area distribution 

Denote links on the lattice pointing right, up, left, and down with operators x, y, x^^, y^^, 
respectively. By convention, all walks begin at the origin. A walk of length A^ is then defined 
by a sequence of links {li, . . . ,l]\f}, where each 4 can be one of the four operators above. For a 
walk to be closed, the number of x's has to be equal to the number of x~^'s, the same for the 
y's and y~^'s, hence an even A^. 

In order to calculate the algebraic area A{li, . . . ,l]\f) enclosed by a walk (positive/negative 
when encircled anticlockwise/clockwise), it is sufficient to note that for two walks differ- 
ing only by an interchange of a pair of subsequent links, A{li, . . . ,lk-i,x,y,lk+2, ■ ■ ■ ^I-n) = 
A{li, . . . , y, X, lk+2, ■ ■ ■ ,In) + 1- Introduce the g-commutator 

xy = qyx , (2) 

where q ^ and q ^ ±1. The full commutation table follows immediately: e.g., multiplying 
both sides by x~^ on the left and on the right yields yx^^ = qx^^y^ etc. Then the algebraic area 
of a closed walk is related to the product of all links constituting that walk as 

= q^(^^'-M _ (3) 

This can be simply understood by defining a normally ordered walk as 

■ ■ ■ ,y'^ ,y,y, ■ ■ ■ ,y,x,x, . . . ,x,x~^ ,x-^ , . . . , x^^} . (4) 

Obviously, the area enclosed by such a walk vanishes, and the product of its links is 1. Calculating 
the area enclosed by an arbitrary walk, then, reduces to normally ordering the product of its links 
and summing up all the powers of q generated by noncommuting links being "carried through" 
each other. E.g., a 1 x 1 square encircled anticlockwise is y~^xyx~^ = qy~^yxx~^ = g, and the 
area is 1. 
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A closed walk of length N containing M instances of x contains the same number of x ^'s 

N 
2 



and, correspondingly, — M of each of y's and y ^'s. The number of such closed walks is 



m 

^^'^ = M!2(f-M)!2 

(A^! permutations of all the links, divided by permutations of identical links), hence the total 
number of closed walks of length is 

= E C^,M = ( . (6) 

We are interested in finding the number Cat (A) of closed walks of length A^ enclosing an algebraic 
area A\ the area probability distribution is then given by 

P.{A) = ^ . (7) 

Oat 

We will search for the generating function of said distribution, 

oo 

ZM= E C^{A)q^. (8) 
yl=-oo 

To elucidate the method, consider at first the set of walks obtained from ^ by shuffling only 
the x's and y's while leaving the x~^'s and y~^'s in place. For a given M as defined above, the 
part of the walk formed by the x's and y's is a "staircase walk" of width M and height 

A^ 

L = --M- (9) 

for brevity, call it an (M, L) staircase walk. Now invoke the g-binomial theorem [ll]: If x and y 
satisfy Eq. jj]), then 

(^+y)^/'= E ^M.L(g)/:c^^ (10) 

M+L=Af/2 

with ZM,LiQ) being the g-binomial coefficient: 



M + L\ ^ [M + L]g] 
L J^- [M],\[L],\ 



ZmA<i)= 7 ] -^^T^, (11) 



where 

L 



[LW. = ll^-^ = lil+q){l+q + q^)---il + q + ... + q'^-'). (12) 
i=l ^ 

The function ZM,L{q) is the generating function of the area distribution of (M, L) staircase walks 
(more precisely, the closures of those walks, obtained by appending x~^ to the end of a walk and 
prepending y~^ at the beginning), in the sense of Eq. ([8|). Indeed, the LHS of (fTOl) is the sum of all 
combinations of products of x's and y's, with any ordering, such that the total number of x's and 
y's is N/2. The number of (M, L) walks within this set that yield a multiplier when normally 
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ordered, i.e., turned into y^x^^ ^ is the coefficient at q"^ in Zm,l{(i)- Respectively, the generating 
function of the area distribution of all staircase walks of length N/2 is Y1iM+l=n/2 Zm,l{q)- 

Consider now rewriting (x + y)^^"^ as {x + y){x + and expand the binomials on both 

sides using Eq. (fTOl) to find 

ZM,Liq) = ZM,L-i{q) + q^ZM-i,L{q) , (13) 

with the initial condition ZQfl{q) = 1. 

There is a simple physical interpretation of ZM,L{q)- Constructing a walk consisting of M 
instances of x and L instances of y amounts to distributing the M x's among the L + 1 slots 
between the y's (slot number being to the right of the last y, slot number A: = 1, . . . , L to the 
left of the fc-th y, counting the latter from right to left). Denote by the number of x's put into 
slot number /c; there is a constraint ^^^=0"^^ ~ normally order the walk, i.e., to move 

all the x's to the right, one has to carry each of the x's through k y's, hence the total area 
is Ylik=o kmk- Now think of the k-th slot as a single-particle state with energy Sk = k and of the 
x's as bosons distributed among those states with occupation numbers rrik. Then the number 
of walks with an algebraic area A is the multiplicity of the multiparticle level with energy A. 
Respectively, the generating function of the area distribution, ZM^iiq)-, is equal to the partition 
function of M bosons in L + 1 single-particle states with energies 0, . . . ,L, with q = e"^/^. The 
recurrence relation (fT3]l can be interpreted as the first addend on the RHS being the sum over 
all multiparticle states in which the highest single-particle level (whose energy is L) is empty; 
the second addend being the sum over all states in which that level is occupied by at least one 
boson. 

Note that levels and bosons can actually be interchanged, since ZM,L{q) = Zi^M{q)- At 5 = 1 
(infinite temperature) relation (fT3l) becomes the Pascal triangle equation, and Zm,l{^) = (^i^) 
is the total number of multiparticle states, or walks. 

Now consider all possible closed walks of length N . In order to obtain the generating function 
of their area distribution, one can generalize ([TU]) as 

(x + y + x~i + y-^ = ZM,,M2M,LMy~'''y'''x''''x-^' ■ (14) 

Mi,M2,Li,L2 
M-i_+M2+Li+L2=N 

Proceeding as above, one concludes that ^M,M,f -M,f -m('?) i« the generating function of the 
area distribution of closed walks containing M instances of x (and as many x ^) and N/2 — M 
instances of y (and as many y~^). For all closed walks of length N , one has to sum over all 
possible values of M: 

N/2 

ZN{q) = ^M,M,f -Af,f -m(9) • (15) 

M=0 

The recurrence relation for ZMi,M2,Li,L2{q) generalizes Eq. (fT3]l : 

ZMi,M2,Li,L2{q) = ZMi,M2,Li-l,L2{q) + ^Mi,M2,Li,L2-l(9) 

+ q^^-^'ZM,-i,M2M,L2iq) + q^''^''ZM„M2~iM,L2iq) > (i6) 

with the initial condition Zq, 0,0, 0(9) = 1- 

The problem of finding the area distribution has been reduced to solving this recurrence 
relation. 
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3 Results 



Obvious symmetry considerations imply Zmi,M2,Li,L2{q) — ■^Li,L2,J\/i,M2(9)) ^is well as 
Zmt,M2Mm(<i) = ^M2,Mi,Li,L2(l/9) (mirror reflection). In general, 



Zmi,M2,Li,L2(.1) = X] Cmi,M2,Li,L2{^)Q'' 



A=-A- 



where the C's are integers and 



A± = max(Mi, M2) max(Li, L2 



I (Ml - M2)(Li - L2)| ± (Ml - M2)(Li - L2) 



(17) 



ZMi,M2,Li,L2i^) is the number of walks involved, i.e., the multinomial coefficient 

(Ml +M2 + Li + L2)! 



^Mi,M2,Li,L2(l) 



Mi!M2!Li!L2! 



(19) 



[Eq. dS]) corresponds to the case Mi = M2 = M and Li = L2 = — M]. At first order in 
9-1, 



7 (\ 7 m\^ (Ml - M2){L, - L2) 
Zhh^WhMMW = Zivh,M2MM\^) ^ 2 - IJ + 

Further, when one of the subscripts vanishes, for example Mi = 0, one hafS 

min(Li,L2) 



(20) 



0,M2,L 



1,^2(9) 



E 

fc=0 



M2 + Li + L2 
k 



M2 + Li + L2 
k-l 



ZM2M-k{<l)ZM2,L2-k{-) ■ 

(21) 

It can be verified by direct calculation, using Eq. (fTSl) . that ^o,M2,Li,L2(9) = ■^o,M2,Li-i,-L2(<?) + 
^o,M2,Li,L2-i(9) + q^^'^'^ Zqm2-im,L2{i)- Furthermore 



(x ^ + y + y^ 



l\N 



-'OM2 



-M2 



(22) 



M2,Lx,L2 
M2+Li+L2=N 



takes place. This generalizes the g-binomial theorem onto the case of three addends. 

Equation (f2T]) relates to the subset of closed walks obtained from (|4| by first moving the 
x's to the end of the walk and then shuffling the x~^'s, y's and y~^'s. Following the same 
line of reasoning as for the staircase walks, the slots between the y's and y~^'s correspond to 
Li + L2 + 1 single-particle states, again numbered from onwards, but now their energies satisfy 
Efc = Ek-i ± 1, where the sign coincides with the sign of the power of y at the k-th position from 
right (the initial condition is still Eq = 0). Hence, the energy of the last state is eLi+L2 = Li — L2- 
The M2 instances of have to be distributed in all possible single-particle spectra stemming 



When q ~ 1, the identity 



E 



M2+L1+ L2 
k 



M2 +L1+L2 
k- 1 



M2 + Li - k\ M2 + L 



M2 



M2 



(M2 + Li +L2)! 
M2\Li\L2\ 



follows. 
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from all possible arrangements of the y's and y~^'s, and Eq. (pTI) is the partition function of M2 
bosons in such a single-particle spectrum, summed over all possible single-particle spectra. 

In the special case Li = L2 = L, the single-particle spectra are in one-to-one correspondence 
with all possible bilateral Dyck paths of length 2L. The partition function in question is then 
related to the statistics of such paths. 

In the absence of a general closed-form solution for ^Mi,M2,Li,L2(<?)) ^Q- ^fT6]l has to be solved 
iteratively (see Table 1). 





iV = 2 


4 


6 


8 


10 


12 


A = 


4 


28 


232 


2156 


21944 


240280 


±1 




4 


72 


1008 


13160 


168780 


±2 






12 


308 


5540 


87192 


±3 








48 


1560 


33628 


±4 








8 


420 


11964 


±5 










80 


3636 


±6 










20 


1200 


±7 












264 


±8 












72 


±9 
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Table 1: Nonzero values of Cn{A) for < 12. 



Note that: 

(i) if Y is even, the maximal possible area is and Cn{^) = N. The maximal area 
is obtained for an anticlockwise square walk x ^), and all such walks are the cyclic 
permutations of {x, . . . ,x,y, . . . ,y, x~^, . . . , x~^,y~^, . . . , y~^}; 

(ii) if Y is odd, the maximal possible area is and Cn{^^) = 2N. The maximal area 
is obtained for anticlockwise rectangular walks {—^ x ^^^) and {^^^ x ^^^j^), and all such 
walks are the N cyclic permutations for one configuration and the cyclic permutations for the 
other one. 

(iii) knowing these values, one can calculate the k-th moment of the area distribution as a 
function of A'' 

00 

Rk{N)= Pn{A)AK (23) 

A=-oo 

The results for R2{N) and Ri{N) coincide with those of Ref. |12| . 

4 Asymptotic limit 

To make the connection with the Harper-Hofstadter model, define a lattice site m = (mi,m2), 
and magnetic translation operators W{m), which satisfy 

W{m)W{m') = W{m + m')e^''^'"i™2-'"2m'i) ^ (24) 

where 7 = 2T^(f)/(j)Q with (f) the flux of the magnetic field per unit cell, the flux quantum. 
Equation ([24l) is the Harper-Hofstadter counterpart of Eq. ([2]) for the algebraic area distribution 
of random walks. 
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The Harper-Hofstadter Hamiltonian is [9l [H] 

H = 



|mi| + lm2|=l 



such that 



TtH 



N 



(25) 



(26) 



where the trace on the LHS is per unit cell and the summation on the RHS is over all closed 
walks of length N. Setting q e'^'/^, the mapping of the problem of random walks onto the 
Harper-Hofstadter problem follows as 



Prom Eqs. one finds that 

Zn{1) 



\'y=x/N 



00 

E 

A=~oo 



PNiA)e 



iAx/N 



(27) 



(28) 



thus establishing, via ([271) . a one-to-one correspondence between the algebraic area probability 
distribution Pn{A) and the A^-th moment of the Harper-Hofstadter spectrum Tv . 



Evaluating Tr H^\^^^/^ in the large N limit one concludes that 



ix/N\ 



4^+^ x/4 
2ttN sinh(x/4) 



1 



2N sinh2(x/4) 



+ 0(l/iV2) 



(29) 



must take place. Obviously, for a finite N, ([29]) cannot hold for all x, since the LHS is periodic 
with period 2TrN. Still, for x < ttN, calculating Z]y from Eqs. (fTSil - lfTBl) results in an excellent 
agreement, even for rather small values of (see Fig. 1). 



5 Conclusion 

The area distribution of closed random walks on a square lattice stems from the noncommutative 
nature of the links that form such walks. This is reflected in the fact that the generating function 
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of that distribution is directly connected with the moments of the spectrum of the Harper- 
Hofstadter Hamiltonian — a sum of four nearest-neighbor magnetic translation operators, which 
become noncommutative in the presence of an external magnetic field. Clearly, the recurrence 
relation ([16]) for the generating function can be interpreted in terms of multi-body partition 
functions, as it has been done for staircase walks. The complexity of the Harper-Hofstadter 
spectrum is encoded in some way in this recurrence relation. 

A generalization of the g-binomial theorem has been obtained in the case of three addends. 
In the case of four addends, closed- form expressions for the generating function might involve 
g-deformed multinomial coefficients — as suggested by the three-addend solution. 
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